Equilibrium and non-equilibrium electron tunneling via discrete quantum states 
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Tunneling is measured via the quantum levels of a metal nanoparticle. We analyze quantitatively 
the resonance energies, widths, and amplitudes, both in the regime where only one state is accessible 
for tunneling and in the non-equilibrium regime when additional states are made accessible one-by- 
one. For tunneling through one state, our results agree with expectations for sequential tunneling, 
but in the non-equilibrium regime the resonances are broadened and shifted in ways that require 
taking into account electron interactions and relaxation. 

PACS numbers: 73.22.-f, 73.23.Hk, 74.80.Bj 



i 

-a 
c 

o 
o 



> 

(N 

o 

O 



i 

C 

o 
o 



In nanometer scale devices, electron transport can oc- 
cur through well-resolved quantum states. Transport in 
semiconductor quantum dots, metal nanoparticles, and 
molecules can all be understood within a similar frame- 
work, in terms of the energies of the states and the rates 
for transitions between states. Here we make a quantita- 
tive investigation of the fundamental processes at work 
when electrons tunnel through a discrete-state system, 
using an Al nanoparticle in a single-electron transistor. 
By varying gate and source-drain voltages (V g ,V), we can 
manipulate electron flow controllably through one state, 
or through many, and we can extract tunneling rates for 
each level. When only one state participates in tunneling, 
the transport properties are in accord with expectations 
for simple sequential tunneling. However, for larger volt- 
ages, the resonance energies, widths, and currents can 
all be modified by the population of excited quantum 
states. For an understanding of the high-I^ regime, non- 
equilibrium transitions beyond those considered previ- 
ously must be taken into account. 

A cross-sectional device schematic is shown in 
Fig. 1(a). The use of an aluminum particle with alu- 
minum oxide tunnel barriers provides mechanical and 
charge stability, and allows V and V g to be varied with- 
out significantly altering barrier resistances. Fabrication 
is done using electron-beam lithography and reactive ion 
etching to create a bowl-shaped hole in a silicon-nitride 
membrane, with a minimum diameter ~ 10 nm. A gate 
electrode is formed by depositing 18.5 nm Al, followed 
by anodization to 3.5 V in an oxygen plasma, and then 
deposition of 8.5 nm of SiO^ |Q. The rest of the device is 
made by depositing a thick Al electrode onto the bowl- 
shaped side of the membrane, oxidizing for 3 min in 50 
mTorr of O2, depositing 1.5 nm Al onto the other side 
of the device to make a layer of Al nanoparticles, oxidiz- 
ing, and then depositing the lower Al electrode. Device 
parameters are determined from the large- ^ structure of 
the 4.2-K Coulomb staircase curve The capacitance 
of nanoparticle to the top electrode is Ct, = 7.9 aF, to 
the bottom electrode Cr = 2.7 aF, and the gate capac- 
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(a) Cross-sectional device schematic. (b) 
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0.06 aF. The sum of the resistances of 



V . A 0.06 Tesla field is applied to drive the Al leads normal, 
(c) / vs. V for different V g in the equilibrium regime. The 
steps correspond to resonances III and IV. Insets show the 
corresponding tunneling transitions. 

the two tunnel junctions is i?s ~ 3Mf2, with individual 
resistances sufficiently large that the intrinsic widths of 
the quantum states are smaller than ksT. Assuming a 
roughly hemispherical |Q particle shape, and a capaci- 
tance per unit area of 0.05 aF/nm 2 fl, we estimate a 
nanoparticle diameter ~ 10 nm. 

In Fig. 1(b), we plot the differential conductance 
dl/dV as a function of V g and V, when the sample is 
cooled in a dilution refrigerator with copper-powder fil- 
ters on the electrical leads. The lines in the figure are due 
to tunneling resonances through discrete quantum states 
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in the nanoparticle. Lines having a positive slope corre- 
spond to tunneling thresholds across the lower-resistance 
junction L, and negative slopes are thresholds across 
junction R jf|. The discontinuity evident in the figure 
is due to a Vg-driven change in the charge on another 
nanoparticle adjacent to the one through which tunnel- 
ing occurs. This merely shifts the electrostatic potential 
of the current-carrying particle. The intrinsic energies, 
current levels, and widths of the resonances are not oth- 
erwise altered, so that the full dl / dV -spectrum can be 
constructed. From the absence of spin-Zeeman splitting 
in a magnetic field (not shown) for resonance lines I and 
II, we can identify these transitions with tunneling from 
an odd number of electrons no on the particle to an even 
number Since these resonances require increased |V| 
as a function of V g , they are uq — > no — 1 transitions. 
(See Fig. 2.) Resonances III and IV correspond to even 
(no — 1) —> odd no transitions. The large gaps in V be- 
tween each of resonances I and II and the next parallel 
lines are due to the energy difference ~ 2A between a 
fully paired superconducting state in the Al particle and 
the next lowest-energy tunneling state with 2 quasipar- 
ticles Q. 

Equilibrium Regime: We first consider the region of V g 
near —650 mV, where V g and V can be adjusted so that a 
single spin-degenerate quantum state (state 0) is accessi- 
ble for tunneling (I-V curves are shown in Fig. 1(c)). For 
the case under consideration, in which the quantum level 
is either empty or singly occupied, the current predicted 
for sequential tunneling |3]-|| is 

~ E (1 + hhaL + (1 + fnhoR 1 ' 

where 7ol (ior) is the bare rate for an electron to tun- 
nel from the quantum state to an unoccupied den- 
sity of states in electrode L (R) and fi is the occupa- 
tion probability for states in electrode i with energy 
equal to the resonance state (for a thermal distribution 
fi = [1 + exp[(eo — Hi)/kT]\ , with eo the energy to 
occupy the quantum state, and /i^ the chemical potential 
in electrode i). Our observations are in excellent accord 
with this model. For instance, the tunneling current 
through the quantum state is not the same for both bias 
directions, being I + = 16.6 ±0.1 pA for V>0 (i.e. cross- 
ing lines II or III) and I- = -8.4 ± 0.1 pA for V < 
(crossing lines I or IV). This has been observed previ- 
ously PJlOll and is a consequence of spin-degeneracy. For 
V > 0, electrons tunnel across the high- resistance rate- 
limiting tunnel barrier into an empty state, so that cither 
spin-up or spin-down electrons can tunnel. For V<0, the 
rate-limiting step is for an electron of a given spin on the 
particle to tunnel through the high-resistance barrier, 
and the current level is approximately cut in half. By 
equating the measured currents to Eq. (1), we determine 
Jor — (5.3±0.1) x 10 7 s _1 for the high-resistance junction 
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FIG. 2. Tunneling diagrams depicting tunneling tran- 
sitions active for resonance lines I-IV in a non-equilibrium 
regime. Black spins represent the ground-state electron con- 
figuration. Black arrows indicate the threshold tunneling 
transition, and gray arrows denote other transitions that con- 
tribute to the current for the value of V depicted. 

and 70L ~ 7 x 10 9 s _1 . In addition, for resonances (e.g., 
Ill, IV) in which the spin-degeneracy of the state is split 
by a magnetic field, the currents through the two Zeeman 
states for a given bias direction are not equal @Jll[], in 
agreement with the simple tunneling theory. The max- 
imum current through the lower-energy Zeeman state is 
e7oi7ofl/(7oL + 7ofl) for either bias direction, and the 
second state then adds current to produce the maximum 
allowed by Eq. (1). 

An interesting feature of Eq. (1) when the tunneling 
threshold is across the lower resistance barrier (peaks II 
and IV) is that the maximum of dl/dV does not occur 
exactly when fa = eo- Instead, the resonance is shifted 
to lower \V\, by an amount oc T, so that it shifts with 
T. This is due to charge accumulation on the particle as 
the Fermi function sweeps by the energy of the quantum 
state, which limits the current at larger |V| on account of 
Coulomb blockade. Our observations of this effect (not 
shown) are equivalent to results of Deshpande et al. [[l0| , 
although the data in JlO[ | were compared to an approxi- 
mation that differed from Eq. (1) [Q. 

Non-equilibrium Regime: We can controllably tune the 
device so that more than one quantum state can partic- 
ipate in tunneling. This is illustrated by following line 
II in Fig. 1(b). This line corresponds to processes which 
are initiated by an electron tunneling off the nanoparti- 
cle from the quantum state to electrode L. However, 
as one follows line II to higher V, past negative-sloping 
resonance lines which intersect line II, these lines indi- 
cate that the subsequent tunneling of an electron from 
electrode R back onto the nanoparticle can proceed via 
many different energy levels other than state 0. The 
total current under these conditions can be modeled 
by a master equation which takes into account all al- 
lowed transitions between the energetically-accessible 
no- and (no — l)-electron states ||. In Fig. 2, we show 
tunneling diagrams depicting representative accessible 
states for resonances I-IV in a non-equilibrium regime. 
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FIG. 3. Properties of the resonance lines as V g is varied 
to change the value of the V at which the resonances appear. 
The second tunneling states become energetically accessible 
for \V\ > 0.18 mV for resonance line I, 0.19 for line II, 0.50 
mV for line III, and 0.64 mV for line IV. (a) Magnitude of tun- 
neling current, (b) Width of the dl /dV peak as a function of 
V, expressed as an effective temperature, (c) Deviation from 
linearity for the ^-position of peak III as a function of V g . A 
line was fit to the peak position in the equilibrium regime be- 
tween V g = —650 mV and —701 mV, and this was subtracted 
from the measured positions. 

In Fig. 3(a) we plot the step height in current associ- 
ated with resonance lines I-IV, as V g and V are tuned 
to follow the lines in the |V|-V^ plane. Peaks I and III 
have approximately constant amplitude, while the cur- 
rents for peaks II and IV grow quickly as |V| enters the 
non-equilibrium regime. This can be understood triv- 
ially. For peaks I and III, the tunneling threshold is 
across the higher-resistance tunnel junction R. This junc- 
tion is always rate-limiting and it matters little how many 
transport channels are available across junction L. For 
peaks II and IV, the tunneling threshold is across the low- 
resistance barrier L, but the rate-limiting process is how 
quickly electrons can tunnel across the other barrier. As 
| V | is increased, more quantum levels contribute to this 
process, and the current grows. By measuring the current 
along peaks II and IV as levels are added one by one, and 
then fitting to master-equation results, we can measure 
rate-limiting tunneling rates for each quantum state: the 
ji R , for i = to 5, are (5.3, 15.7, 8.0, 16, 15, 9) x lO 7 * -1 , 
7i « 3 x 10 9 s _1 , and (assuming relaxation effects are 
negligible) 7_ir s» 17 X 10 7 s _1 . 

Figure 3(b) shows the widths of the tunneling reso- 
nance lines I-IV. These were determined by fitting the 



V dependence of each conductance peak to the deriva- 
tive of a Fermi function, and then converting the voltage 
width to effective temperature by multiplying by the ca- 
pacitance ratio (e/fcs)CL/(CL + Cr) for peaks I and III 
or (e/k B )C R /(C L + C R ) for peaks II and IV g). In either 
the equilibrium or non-equilibrium regimes, the predic- 
tion of the simplest master equation |^] is that the peak 
shape should be a derivative of the Fermi function with 
a width approximately equal to the electronic temper- 
ature in the electrodes. Our measurements agree with 
this model within the equilibrium regime, with a con- 
stant electron temperature T w 90 mK. This is signifi- 
cantly higher than the T achieved in non-gated tunneling 
devices using our apparatus, 45 mK, and we ascribe the 
difference to heating by leakage current from the gate to 
electrode R. Line III is broader than the others at low 

V because the magnetic field of 60 mT applied to drive 
the Al electrodes normal produces an unresolved Zeeman 
splitting (AE/k B = 2fi B H/k B = 80 mK). 

As |V| is increased into the non-equilibrium regime, 
peaks I and III undergo large increases in width and peak 
IV broadens slightly, while peak II shows no measurable 
change. The differences are not merely an effect of heat- 
ing in the electrodes, because peaks II and IV have the 
largest magnitudes of current and power. We suggest 
that these measurements can be explained as a conse- 
quence of electronic interactions in the non-equilibrium 
regime, by a mechanism due to Agam et al. ]l2] ]. Con- 
sider resonance line III, for which the tunneling thresh- 
old corresponds to an electron entering quantum state 
0. For V > 0.50 mV the next tunneling event, which 
discharges the particle, may occur out of different, lower 
energy states (see Fig. 2), leaving an electron- hole exci- 
tation on the nanoparticle. Agam et al. suggested that if 
this non-equilibrium state does not relax before the next 
electron tunnels onto the particle, it can shift the energy 
of the tunneling resonances on account of an alteration 
of the electron-electron interaction energy In past 

work on smaller aluminum particles, shifted transitions 
were resolved individually [jl2"f ; however the relative shift 
is expected to decrease with increasing nanoparticle size 
p2| , so it is reasonable that the shifts would produce 
only broadened resonances for the 10 nm particle un- 
der investigation here. Because a growing ensemble of 
different non-equilibrium states can be excited with in- 
creasing |V|, this mechanism can explain the increase in 
width of line III as a function of \V\. The same non- 
equilibrium mechanisms should also come into play for 
line IV, for V < —0.64 mV, but the broadening here is re- 
duced because the threshold tunneling event is across the 
lower-resistance junction, L. Barrier R quickly becomes 
rate limiting as line IV is crossed, so that higher-energy 
non-equilibrium resonances do not add significant addi- 
tional current. In order for the non-equilibrium mech- 
anism to apply for line III, the relaxation rate of some 
non-equilibrium excitations to the ground state must be 
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comparable to or slower than jqr — 5.3 x 10 7 s~ 1 . The 
rate predicted by Agam et al. for spin-preserving energy 
relaxation in aluminum particles is ~ 10 8 s _1 |l2[ . 

Resonances I and II are a different case, because 
the tunneling threshold corresponds to an electron leav- 
ing quantum state 0. The subsequent tunneling event, 
adding an electron back to the nanoparticle, may for large 
V occur into higher-energy states (see Fig. 2(a,b)), but 
nevertheless this excitation alone cannot produce a non- 
equilibrium shift in the energy of subsequent discharg- 
ing transitions. The reason is that only this electron is 
free to tunnel off the nanoparticle; there is no electron 
in quantum state whose transition energy might be 
shifted. Therefore within the picture of Agam et al. [fl2|| , 
no non-equilibrium broadening should be expected for 
levels I and II, in conflict with the data for level I. This 
discrepancy can be explained if non-equilibrium excita- 
tions on the nanoparticle can be generated not only by 
the tunneling transitions on or off the particle that have 
been considered previously , but also by transitions in 
which a high-energy electron relaxes within the nanopar- 
ticle and produces an electron-hole excitation. Broaden- 
ing would then be generated by the Agam mechanism. 
Within this scenario, the difference between the broad- 
ening visible for resonance line I and the lack of broad- 
ening of line II would follow from the fact that for peak 
II a high-energy electron on the particle can quickly exit 
through the low-resistance tunnel junction L, while for 
peak I the high-energy particle must exit through the 
high-resistance junction R, giving a much longer resi- 
dence time during which relaxation transitions can occur. 
In order for line I to be broadened, the fastest relaxation 
rates must become comparable to jqr = 5.3 x 10 s _1 as 
\V\ increases. 

By tuning V g and V into the non-equilibrium regime, 
the apparent energies of the dl/dV peaks can also be 
changed. This is clearest for line III (Fig. 3(c)), which un- 
dergoes a shift of 33/iV to lower voltage when the thresh- 
old for non-equilibrium tunneling via state -1 (line V in 
Fig. 1(b)) is crossed. Because we have measured rate- 
limiting tunneling rates for the energetically-accessible 
states from the current amplitudes, we can test whether 
this shift can be explained by the simplest master equa- 
tion M, which assumes that the underlying energies of 
the quantum states are not changed by non-equilibrium 
interactions. Only one relevant parameter is not deter- 
mined previously: x — 7ol/7-il- The solution of the 
master equation does predict a voltage shift (oc T) for 
the conductance peak compared to the equilibrium case, 
and for x < 0.15 it can explain the full value of the 
experimental shift. However, we judge this to be im- 
probable, because the measured values of "fm fall within 
a more narrow distribution. For x ~ 1 in Eq. (2), the 
predicted shift is 15/iV - much smaller than we mea- 
sure. We can more naturally explain the full value of 
the shift by again taking into account that the presence 



of a non-equilibrium excitation can change the energy of 
a tunneling transition. In the equilibrium regime, the 
occupation of quantum state corresponds to a tran- 
sition from a fully-paired superconducting state on the 
aluminum particle to a state with one high-energy quasi- 
particle; in the non-equilibrium case, the transition can 
be from a state with two quasiparticles to one, with a 
transition energy lowered by ~2A«0.35 meV [Q. This 
is much bigger than fc^T ~ 10/ieV, and in this case the 
observed resonance is shifted to lower \V\ by an amount 
oc T because electrons in the tail of the Fermi distribution 
can excite the non-equilibrium state and open the lower- 
energy current channel. For x — 1 and T ~ 90 mK the 
master-equation result is that the measured shift can be 
produced by a non-equilibrium lowering of the transition 
energy by any amount greater than 20/ieV^2kBT. 

In summary, we have made a systematic study of the 
transition energies, widths, and current levels for elec- 
trons tunneling through discrete quantum states in a 
nanometer-scale electronic device. When transport oc- 
curs through a single quantum level, the results are in 
agreement with the expectations of sequential tunneling. 
At large values of |V|, the non-equilibrium population of 
excited electronic states, together with electron-electron 
interactions, can modify the widths and apparent ener- 
gies of the tunneling resonances. 
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